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Abstract 

In this study, we consider the asymptotic behaviour of the first discrete Painleve equa¬ 
tion in the limit as the independent variable becomes large. Using an asymptotic series 
expansion, we identify two types of solutions which are pole-free within some sector of the 
complex plane containing the positive real axis. Using exponential asymptotic techniques, 
we determine the Stokes Phenomena effects present within these solutions, and hence the 
regions in which the asymptotic series expression is valid. From a careful analysis of the 
switching behaviour across Stokes lines, we find that the first type of solution is uniquely 
defined, while the second type contains two free parameters, and that the region of validity 
may be extended for appropriate choice of these parameters. 


1 Introduction 


We consider the asymptotic limit n —> oo of the first discrete Painleve equation (dPi) 


an + /3 

w n+1 + w n + w n -i = -b 7 , 

w n 


( 1 . 1 ) 


and apply methods of exponential asymptotics to describe solutions which demonstrate asymp¬ 
totically pole-free behaviour within certain sectors. This nonlinear equation also appears in the 
study of orthogonal polynomials, known there as the Freud equation mils]. Previous studies 
mm considered the far-field behaviour of solutions to dPi, finding the asymptotic solution in 
terms of elliptic functions. Joshi 0 and Vereshchagin [6; also determined discrete analogues of 
asymptotically pole-free solutions behaviours by Boutroux for the continuous first Painleve equa¬ 
tion (Pi). In this paper, we find the full formal series solutions including exponentially-small 
terms, study the influence of Stokes phenomena on these solutions, and deduce their regions 
of validity. The resultant asymptotic solutions share several features with the asymptotically 
pole-free solutions to Pi known as tronquee and tritronquee solutions. 

Equation (11.11) tends to Pi in the scaling limit w n = 1 — h(2y{x) — c/2), x = \fhn — 6/2 — c 2 /24, 
with a = h 3//2 a , /3 = —7 + 3 + h 2 b, 7 = 6 + he, as h —> 0. Pi is one of the six Painleve equations 
which play a key role in a range of physical problems, in particular arising as universal models 
in reductions of the nonlinear Schrodinger equation [ 7 ], the Korteweg-de Vries equation mm 
and the Boussinesq equation TO], and appear as models in such diverse fields such as statistical 
mechanics nn, nonlinear optics [ 12 ] , electrochemistry m and plasma physics [5:. The first 
discrete Painleve equation also arises in physical applications, such as the study of quantum 

gravity [HUS]. 

Motivated by these applications, the study of the Painleve equations in asymptotic limits has 
been widely pursued. However, the corresponding study of discrete Painleve equations remains 
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sparse. In this paper, we apply techniques of exponential asymptotics for differential-difference 
and difference equations developed by King and Chapman |16] in their analysis of the Frenkel- 
Kontorova model for atomic lattices to dPi, in order to obtain asymptotic approximations which 
are accurate up to exponentially-small correction terms. 

We note that there have been some previous investigations into exponential asymptotics for 
difference equations, such as that of [m m- The study by Olde Daalhuis na performed a 
Borel-plane analysis to consider the behaviour of second-order linear difference equations and 
to determine exponentially-improved solutions in terms of an inverse factorial series. In fact, 
these Borel-plane methods could be applied to the current investigation, and controlled error 
estimates for the asymptotics could be obtained in this fashion. The primary reason for using 
the Stokes-smoothing technique described in m is that the analysis may be performed directly 
on the nonlinear difference equation, rather than requiring the solution to be formulated in terms 
of an integral expression; this feature permits the analysis to be easily extended to a wide range 
of more complicated problems. 


1.1 Background 


The asymptotic behaviour of solutions to the Painleve equations was studied by Boutroux m, 
who considered the behaviour of the first Painleve equation (Pi), 


d 2 y 
dx 2 


6 y 2 + x, 


in the limit |x| —» oo. Boutroux identified special solutions which are asymptotically free of poles 
in particular sectors of the complex plane of angular width 47r/5, and furthermore discovered 
five unique solutions that are pole-free in sectors of angular width 87r/5. These solutions are 
known as tronquee and tritronquee solutions respectively, and tend tow~ ± \J—x/d as |x| —> oo 
within the asymptotically pole-free sectors. In [20], the real tritronquee solution was shown to 
have no poles on the whole semi-axis x < 0. Several asymptotic studies of Pi have been carried 
out in These asymptotic approximations have played essential roles in many of the 

physical problems described earlier mm- 


1.2 Discrete Painleve equations 


In this study, we are interested in the asymptotic behaviour of solutions of the first discrete 
Painleve equation. The discrete Painleve equations are a set of integrable nonlinear difference 
equations, which tend to the ordinary Painleve equations when the continuum limit is taken. 

The discrete Painleve equations first appeared in the study of solutions to the continuous 
Painleve equations. It was discovered that these equations possess Backlund transformations, 
which map one solution to another solution with different parameter values, or even to a solu¬ 
tion of a different Painleve equation (see for example [TO] 115)- These transformations have been 
used to construct hierarchies of solutions to the Painleve equations, often in terms of rational 
or special functions. 

It was noted by Fokas, Its and Kitaev [14; [15] that dPi arises in the study of Backlund transfor¬ 
mations associated with the fourth Painleve equation (Piv)- This equation has two parameters, 
a and 6, and is given by 


d 2 w 

cb 7 


/\2 


K) 

2 w 


3 w 2 , o o , b 

-b 4 zw 2 + 2 (z 2 - a)w + —. 

2 w 


If we associate a choice of parameters ( a ni b n ) of Piv with a solution denoted by w n (z), a 
hierarchy of solutions may be found corresponding to parameter choices 


a n = -§Ai(-l) n + \J—2b n = //(-!)" + n, 


( 1 . 2 ) 
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where p is an arbitrary complex-valued constant. The recurrence relation for w n is subsequently 
given by 

(Wn+l + Wn + W n - i)w n = -2 ZW n - U + fl, (1.3) 

which is simply a version of dPi with appropriately chosen parameters. A variant of this re¬ 
lationship was used by [23] to describe hierarchies of solutions to Piy in terms of parabolic 
cylinder functions, as well as in uni to construct solutions in terms of rational functions, as well 
as Okamoto and Hermite polynomials. By studying the asymptotic behaviour of the discrete 
Painleve equations, we are able to better understand the behaviour of these solution families. 

The first discrete Painleve equation has also featured in the study of problems arising outside of 
the field of integrable systems, particularly relating to the study of orthogonal and orthonormal 
polynomials 12113- A set of orthonormal polynomials p n (x), where n > 0, is defined with respect 
to some weight function w(x) as satisfying 

/ OO 

Pn(x)p m (x)w(x)dx = S m>n , (1.4) 

-OO 

where S mi „ = 1 if m = n and is equal to zero otherwise. In general, these polynomials satisfy a 
recurrence relation 

xp n {x) = a n+1 p n+ i(x) + b n p n (x) + a n p n - i(x), 

where a n and b n themselves satisfy particular recurrence relations, described for example in [2|. 
In the case of orthonormal polynomials with respect to the Freud weight, ui(x ) = exp(— kx 4 — 
px 2 ), we find that b n = 0 (this holds for all even weight functions), and that u n = a\ satisfies 

4:Ku n (u n -i +u n + u n+ 1 ) + 2pu n - n = 0, (1.5) 

which, again, is dPi, with appropriately chosen parameters. By analyzing the solutions to dPi, 
we can determine the behaviour of families of orthogonal polynomials, even in the limit that 
their degree becomes large. It is this relationship to orthonormal polynomials that underpins the 
appearance of dPi in the study of quantum gravity models mm- We also note in passing that 
the variation of w n with respect to the parameters n and p can be computed directly from the 
fourth continuous Painleve equation with appropriate choice of parameters [2], thus confirming 
the existence of another connection between dPi and Piv- 

In recent decades, however, discrete integrable systems such as dPi have become an object 
of study in themselves. In part, this is due to these equations being convenient examples of 
nonlinear difference equations which, due to their integrability, are more tractable than many 
equations within this class. Hence, they provide a useful testing ground for mathematical 
techniques intended for wider applicability. 

The primary motivation for studying these equations, however, is to better understand their 
solutions. These equations are important discretizations of the continuous Painleve equations; 
their particular significance is due to the fact that they retain the integrability property upon 
discretization. They arise as independent entities in applications and possess interesting solu¬ 
tions, whose properties permit us to understand the behaviour of classes of functions, such as 
the orthonormal polynomials described above. While addressing such questions about all solu¬ 
tions lies beyond the scope of the present study, it motivates our interest in understanding the 
asymptotic behaviour of the discrete Painleve equations. 

1.3 Exponential asymptotics and Stokes lines 

The formal series we study in this paper are shown to be divergent (see Section [2]) and contain 
exponentially-small terms that are invisible to conventional asymptotic approaches. The central 
idea of exponential asymptotic methods is that a divergent asymptotic series, when truncated 
optimally, produces an error term that is exponentially small in the asymptotic limit [26 1. The 
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problem may then be rescaled to directly determine the behaviour of this remainder term, 
allowing the exponentially-small component of the solution to be determined separately to the 
asymptotic series itself. This idea was introduced by Berry [27 28. 29 , and Berry and Howls 
EDI. who used these methods to determine the behaviour of special functions such as the Airy 
function. 

As we uncover this exponentially-small behaviour, we will discover two varieties of important 
curves, identified originally by Stokes m and known as Stokes lines and anti-Stokes lines. 
These curves relate to the behaviour of exponentially-small components of the solutions. As 
Stokes lines are crossed, the exponentially-small behaviour experiences a smooth, rapid change 
in value in the neighbourhood of the curve itself; in many cases, the behaviour may switch 
off entirely on one side of the Stokes line. Hence, Stokes lines allow us to determine which 
exponentially-small contributions may be present in certain regions of the plane. Anti-Stokes 
lines determine directions where exponentially-small behaviour changes to become exponentially 
large, and therefore mark the boundaries of the region in which the asymptotic series expression 
for the solution is valid. 


In this investigation, we will apply an exponential asymptotic method developed by Olde Daal- 
huis et al. [32) for linear differential equations, and extended by Chapman et al. [33] for 
application to nonlinear differential equations. We provide a brief outline of the process here, 
however more detailed explanation of the methodology may be found in these studies. 


To optimally truncate an asymptotic series, we require a general form for the asymptotic series 
coefficients, which is often algebraically intractable. In singular perturbation problems, Dingle 
[34] noted that successive terms in the asymptotic series expansion are obtained by repeated 
differentiation of an earlier term in the series. Hence, the late-order terms a m of the resultant 
asymptotic series typically diverge as the ratio between a factorial and the increasing power of 
a function y, such as 


Q"m 


AT(m + 7 ) 

x m +1 


as m —>■ 00 , 


( 1 . 6 ) 


where T is the gamma function defined in [35] , and A, 7 and y are functions that do not depend 
on to. In particular, the singulant y is equal to zero when the leading-order behaviour is singular, 
ensuring that this singularity is present in all higher-order terms. Chapman et al. [ 33] therefore 
propose a factorial-over-power ansatz m for the late-order terms, which may then be used to 
optimally truncate the asymptotic expansion. 


The next step, described by Olde Daalhuis et al. [33], is to substitute the truncated series 
expression back into the original problem, obtaining an equation for the exponentially-small 
remainder term after truncation. The positions of Stokes lines are found by determining the 
curves along which the remainder term varies rapidly, and the behaviour of this remainder is 
found by solving the remainder equation in the neighbourhood of these Stokes lines. 


Away from the Stokes line, the governing behavior is obtained by solving the homogenous 
version of the remainder equations, generally giving an expression of the form Ae~ x ^ e (this is a 
consequence of the analysis in Chapman et al. [33], although it is also shown as a part of the 
Stokes smoothing analysis contained herein). Stokes switching will occur when the switching 
exponential is maximally subdominant; this occurs when the singulant is purely real and positive. 
Anti-Stokes lines may be found when the exponential term switches from being small to large; 
this occurs when the singulant is purely imaginary. These are therefore convenient conditions 
to determine the possible location of Stokes and anti-Stokes lines, as observed by Dingle 34] . 
and we apply them to the series solutions of dPi. 


1.4 Paper outline 

In Section 2, we find the full formal series expansions of the (pole-free) solutions of dPi and 
provide recurrence relations for the coefficients. It is important to note that we find two types of 
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series solutions, called Type A and B solutions. We also describe the late-order terms in these 
series. We determine the Stokes structures of the Type A solutions in Section 3, which involves 
descriptions of the exponentially small terms in the solutions. These are used to find the regions 
in which the solutions are pole-free. We carry out the analysis of Type B solutions in Section 4. 
Finally, we discuss the results and conclusions of the paper in Section 5. 


2 Asymptotic series expansions 


In this section, we expand the solution as a power series in the limit n — > oo, and determine a 
recurrence relation for the coefficients of the series. Finally, we determine a general expression 
for the late-order terms of the series. 

Consider the first discrete Painleve equation: 

an + P 

W n + 1 + Wn + W n - 1 = -- + 7 . 

Wn 

Now, simply scaling into the far field requires the behaviour to be smooth in this limit, however 
we also wish to consider asymptotically pole-free solutions that allow for the odd and even terms 
to vary rapidly. Hence, we consider separately the odd and even terms in the series, such that 


w n 


y n if n even 
z n if n odd ’ 


( 2 . 1 ) 


where y n and z n do not necessarily have the same asymptotic behaviour as n — > oo. Equation 
(I2.2al) may now be represented as a system, giving 

(z n+1 + y n + z n -i}y n = an + ft +'yy n , n even (2.2a) 

{y n +i + z n + y n -i)z n = an + p +'yzn, n odd. (2.2b) 


In subsequent analysis, we assume n £ C with iteration n H > n+ 1 acting on domains, and takes 
an initial value as given by a function f(n) in a domain that includes the initial point. This is a 
standard interpretation [36] J37], and underlies the theory of difference calculus. This allows us 
to analytically continue the solution in the complex n-plane. 

We now apply the method of multiple scales in terms of the fast variable n, and a slow variable 
s = en. We scale y n and z n and express them such that 


_ u{s,e) _ v{s,e) 

y™ e i/2 ’ Zn e i/2 > 

which gives the rescaled equations 

(v(s + e, e) + u(s, e) + v(s — e, e))u(s, e) = as + ep + e 1 / 2 q u(s, e), 
(u(s + e, e) + v(s, e) + u(s — e, e))v(s, e) = as + ep + e 1 ^ 2 "fv(s, e). 


(2.3) 


(2.4a) 

(2.4b) 


We assume that u(s,e ) and v(s,e) are analytic in the variable s, and expand the solutions to 
(12.4a.D - (12.4bl) as a Taylor series in s to give the system 


*0,e) +Y^ 


j =o 


^(s,e) + J2 


3=0 


2 e 2j u( 2 ^) (s, e) 

m 

2 e 2 - 7 rd 2j ) (s, e) 

m 


i(s, e) = as + ep + e 1 ^ 2 7 u(s, e), 
;(s, e) = as + ep + e 1 ^ 2 7 v(s, e). 


(2.5a) 

(2.5b) 
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( 2 . 6 ) 


We now expand the solutions u(s, e) and v(s, e) as a formal power series in e 1 / 2 , such that 

oo oo 

u(s,e) ~ ^2 e m/2 Mm(s), v(s, e) ~ ^ e m/2 v m (s). 


m—0 


m —0 


Equations (I2.5al) - (l2.5bl) therefore become 


OO OO 


OO ^ On oo oo 


EE e ( m + r )/ 2 umur _l_ 


ESEE e^)^^) = as + e(3 + 7 £ e^/ 2 

m=0 


m—0 r=0 
oo oo 


(27)! 

j'=0 v ^ ' m—0 r—0 

oo 2 2j 00 00 


^m(^) 5 


E E e (m+r)/2 ^mWr. + ^ y^-7-y ^2 E e(m+r,)/2 ^m4 2,7) = as + e/3 + 7 E e (m+1)/2 v m (s). 

j=0 ^ J ' m=0 r=0 


m—0 r=0 


m—0 


Matching orders of e as e —> 0 gives 


O(e 0 ) : 

u 0 (u 0 + 2v 0 ) = as, 


Vo(v 0 + 2 uo) = as, 

0(e 1/2 ) : 

2(u 0 ui + u 0 vi + uiv 0 ) = 7 u 0 , 


2(v 0 v 1 + VQUi +viu 0 ) = ryvo, 

Oie 1 ) : 

u\ + 2 u 0 u 2 + 2uqv 2 + 2vqu 2 + 2u\V\ = /3 + 'fui 


v\ + 2v 0 v 2 + 2vqu 2 + 2u 0 v 2 + 2viui = /? + 717 . 


This gives two sets of solutions. The first of these is given by 
uq = —vq = ±\/— as, 
while the second is 


7 4/3 — 7 2 

u 1 =v 1 = u 2 = -v 2 = T 0 ,- 

2 8a/— as 


u 0 = Vo = ± 


u l = Vi = 


1 

6’ 


U 2 = V 2 = ± 


(12/3 + 7 2 )v / 3 
72Vas ’ 


(2.9) 


( 2 . 10 ) 


where the upper and lower choice of signs correspond. In subsequent analysis, we will refer to 
solutions with the leading-order behaviour given by (GU) as Type A, while those with leading- 
order behaviour given by (12.101) as Type B. 

We see that in both solution types, u o and uo contain a singularity at s = 0, and that this 
singularity increases in strength by one half in each subsequent term of the asymptotic series. 
Consequently, we expect that the series will diverge, and exponentially-small Stokes switching 
behaviour will be present in the solution. 

In general, we find 


2m 

[m/2J 


2m—4 j 



3=0 

-i+E 

3=0 

Z 

m 

E 

r—0 

(23) 

U r V 2 m - r -4j = 7«2m-l, 

(2.11a) 

2m 

Lm/2j 


2m—4 j 



E V 3 v 2m- 

i=o 

-i+E 

3=0 

z 

m 

E 

r=0 

(23) 

V r u 2m-r-4 j = 

(2.11b) 


In order to determine the optimal truncation point, we must determine the form of the terms 
in the asymptotic series as n —> oo. We therefore follow .33], and define a late-order ansatz in 
the limit n —> oo: 


U{s)T{m/2 + k) V{s)T(m/2 + k) 

y ( s'j rn /2+k ’ Vm ^(g'jm/2+k ’ 


( 2 . 12 ) 
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where x is the singulant, satisfying x(0) = 0, U and V are the prefactors, and A: is a constant. 
Applying this ansatz to (l2.llaM2.Hbll and matching leading orders of m as m —> oo gives 



(2.13b) 


(2.13a) 


Continuing to the subsequent order as m —> oo, we find u\ = v\ = 7 /2 and u\ = 17 = 7/6 for 
Type A and B solutions respectively, which are always satisfied. At the next order, we find 



(2.14a) 



(2.14b) 


To determine the behaviour of the prefactor and singulant terms, we must consider Type A and 


B solutions separately, and solve (12.13aD - (l2.13bl) to find the singulant x(s), and then (I2.14afi - 


(12.14bl) to find the prefactors U(s) and V(s). 


3 Stokes structure of Type A solutions 


In this section, we solve the late-order term equations for Type A solutions. We use the resultant 


expressions to determine the Stokes structure, as well as the form of the exponentially-small 
components of the solution. From these results, we determine the region of the complex plane 
which is asymptotically pole-free. 

3.1 Late-order terms 

The singulant equations given in (l2.13aM2.13bl) reduce to 



(3.1) 


Recalling that we are considering the late-order terms in the limit m —> 00 , the leading-order 
behaviour of x may be found by taking the series to be infinite, introducing error into the 
singulant behaviour that is only exponentially-small in this limit. This gives 


cosh(— x') = 0 . 


(3.2) 


The solutions to this expression take the form 


X = ±— +M 7ri, 


(3.3) 
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where M is any integer. In general, the late-order behaviour of u n and v n will be the sum of 
terms associated with the singulants given in m, with each term associated with a particular 
choice of sign and M. However, in a similar fashion to the analysis carried out in m , we note 
that this sum will be dominated in the n —> oo limit by the two terms associated with M = 0. 
Hence, we may restrict our attention to these components in our analysis. 


Recalling that uq and vo are singular at s = 0, and hence that x(0) = 0, we find that the relevant 
singulants are given by 


17TS 17T S 

xi - —! X2 - Y' 


(3.4) 


To find the prefactor associated with these singulants, we solve (I2.14al) - (l2.14bl) . We take the 
limit of the summation terms as n —> oo. Equation (I2.14bll becomes 


2{v 2 V + u 2 V + v 2 U) + 2vo U' sinh(-x') - ^ + 2 v 2 U (cosh(-x') - 1) = 0. 

X 

Applying the form of the singulant given in (13.311 . we find that this reduces to U' = 0. Hence we 
set the prefactor U ( s) = A, where A takes constant value. Recalling that there are two distinct 
singulant contributions, we conclude that there are two constants associated with U, denoted 
as Ai and A 2 . A similar analysis performed on (12.14al) produces V' = 0. Hence, we find that 
there are also two constants associated with V , which we will denote as A 3 and A4. 

Finally, we note that for the ansatz (12.1211 to be consistent with the leading-order behaviour 
given in (EH), we must have k = —1/2. This gives a singularity of the appropriate strength in 
the limit n —>■ 0. Hence, the late-order terms are given by 


AiT(to/2-1/2) , A 2 T(m/2 — 1/2) A 3 T(m/2-l/2) i A 4 T(to/2-1/2) 

Um ~ (i7T S /2) m /2-l/2 + (-i7T S /2) m /2-l/2’ Vm ~ (i7T S /2) m /2-l/2 + (-ms/2) m / 2 ~ 1 / 2 ' 

(3.5) 

The only component of the late-order terms that remains to be found is the value of the constants 
A^ which may be determined numerically. We illustrate this on an example problem in Appendix 

EH 

3.2 Stokes structure 

Having calculated Xi and \ 2 , we recall from Section 11.31 that the behaviour of the exponential 
term will be proportional to exp(—x/e); this is also demonstrated in the full Stokes switching 
analysis contained within Appendix [B] This term will be exponentially small when Re(x) > 0, 
and exponentially large when Re(x) < 0. We also recall that Stokes lines exist where Im(x) = 0 
and Re(x) > 0, while anti-Stokes lines follow Re(x) = 0. Figure HOI illustrates the behaviour of 
the singulants in each of the four quadrants within the complex s-plane. 

We see that there is an anti-Stokes line following the positive real axis, while Stokes lines follow 
the negative and positive imaginary axis for xi and X 2 respectively. To determine the behaviour 
that occurs as the Stokes lines are crossed, we note that the remainder term is of algebraic order 
along anti-Stokes lines. However, this would make it comparable in size to the leading-order 
behaviour on the positive real axis, which is not possible. Hence, we conclude that the remainder 
terms associated with xi and X 2 are not present along the real axis at all, and hence cannot be 
active in the first or fourth quadrant. 

Consequently, we conclude that the exponential contributions associated with xi and X 2 are 
present on the left-hand side of their Stokes lines, but are switched off as they cross from the left 
to the right-hand side (into the fourth and first quadrant respectively). Figure [3~2l illustrates 
the regions in which the exponential contributions are active. 



















Im(s) 



Figure 3.1: Behaviour of the singulants for Type A solutions in each of the four quadrants. The dashed curve 
indicates where the real part of both singulants is zero (and hence corresponds to an anti-Stokes line), 
and the solid curve indicates where the imaginary parts are zero, corresponding to potential Stokes 
lines. The black jagged line represents a branch cut. Recalling that Stokes switching may only occur 
when Im(x) = 0 and Re(x) > 0, we see that xi switches across the curve arg(s) = — 7r/2, while X 2 
switches across the curve arg(s) = n/2. 


Im(s) 

® t 


Xi: Inactive 

X 2 : Active (exp. small) 


Xi: Inactive 
X 2 : Inactive 


Stokes line - > 

- 

Stokes line -► 


Anti-Stokes line 
------—> Re(s) 


Xi: Active (exp. small) 
X 2 : Inactive 


Xi: Inactive 
X 2 : Inactive 


©I 

(a) Regions containing exponential behaviour 


Im(s) 



Figure 3.2: Figure (a) illustrates the regions of the s-plane in which exponential contributions associated with 
Xi and X 2 are active. Switching behaviour of exponential contributions occurs across the Stokes 
lines, denoted as solid black lines, where the associated singulant is indicated by the nearby circled 
numbers. The anti-Stokes line is represented as a dashed curve. The exponential contributions cannot 
be present along the anti-Stokes line, as they would not be exponentially small here, and hence would 
be visible to the algebraic power series. Figure (b) shows the exponentially-small remainder terms 
present in each region. 
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Importantly, we see that as the exponential contributions associated with \i and X ‘2 are present 
only in the third and second quadrants respectively, in which they are exponentially small (as 
Re(x) > 0), they do not cause any difficulty for the validity of the asymptotic expansion. Hence, 
the solutions satisfying (12.91) are valid on the entire complex plane, with the exception of the 
branch cut. We therefore conclude that the solutions are asymptotically pole-free in the entire 
far-held, aside from the negative real axis. 

Finally, now that we know the nature of the switching behaviour as the Stokes lines are crossed, 
we are in a position to calculate the precise form of the remainder term associated with the late- 
order terms, and hence determine the exponentially-small contribution present in the problem. 


3.3 Stokes smoothing 

To determine the behaviour in the neighbourhood of the Stokes line, we return to (I2.5all - (l2.5bll . 
Now, we express u and v as asymptotic expansions truncated after N terms, with remainder 
Rn and Sn respectively, such that 

N N 

u(s,e) = ^2 e m/2 u m (s) + R N {s), v(s, e) = ^ e m/2 v m (s) + S N (s). (3.6) 

m=0 m—0 

If N is chosen so that the series are truncated after their smallest terms, the remainder terms 
will be exponentially small. By considering the behaviour of these remainder terms, we can 
obtain the exponentially-small corrections to the asymptotic series, and use these to determine 
the region in which the asymptotic series provides a valid approximation of the exact solution. 

A full analysis of the remainder terms Rn and Sn is given in Appendix [B] This analysis is quite 
technical; we therefore summarize the results in the remainder of this section. 

From Appendix m we find that the optimal truncation point is given by N opt = 2|y|/e + w, 
where 0 < to < 1 is chosen so that N is integer-valued. The remainder terms given in ESI)-dESl) 
are shown to take the form 

Rn ~ SiAi e _i7rs / 2e + S 2 A 2 e ins ' 2e = <SiA, (-i) s / £ + S 2 A 2 (i) s / £ , (3.7) 

S N ~ SiA 3 e~ ins / 2e + S 2 A 2 e ins ' 2e = SiA 3 (-i) s/e + S 2 A 4 (i) s / e . (3.8) 

where Si is the Stokes-switching parameter that varies rapidly in the neighbourhood of Stokes 
lines. From the analysis in Appendix [Bj we find that the value of Si changes smoothly from 
Si = 4 to Si =0 when crossing from the third to the fourth quadrant, and the value of S 2 changes 
smoothly from S 2 = 0 to S 2 = 4 when crossing from the second to the first quadrant. Hence, 
the exponential behaviour of the asymptotic approximation has been completely determined, 
and is shown in Figure l3T2f b'). 

Consequently, the optimally-truncated asymptotic series for u(s, e) and v(s, e) takes the form 

N Q pt 

u(s , e) ~ ^ e m/2 « m (s) + SrAi (-i) s / £ + S 2 A 2 (i) s / e , 

772 = 0 
N opt 

v(s, t m/2 Vm{s) + S 4 A 3 (-i) 5 / £ + S 2 A 4 (i) s /% 

772 = 0 

where the leading orders of u n and v n are give in & and the late-order behaviour is given in 
m- This expression is therefore an accurate asymptotic approximation up to exponentially- 
small terms, valid in the entire complex s-plane. 


(3.9) 

(3.10) 
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4 Stokes structure of Type B solutions 


In this section, we determine the late-order asymptotic terms associated with Type B. These are 
then used to determine the exponentially-small contributions to the asymptotic solution, and 
subsequently the Stokes structure. From these results, we determine the region of validity for the 
general solution, and show that for particular choices of the two Stokes-switching parameters, 
the region of validity may be extended. 


4.1 Late-order terms 


When considering Type B solutions, we find that the asymptotic series expansion for u and v 
are identical, and hence that u m = v rn for all values of m. We therefore set V(s) = U(s). 

Now, using (12.101) . both singulant equations (I2.1dal) - (l2.1dbl) reduce to 

2 + cosh(— x') = 0, (4.1) 


where we have again taken the upper bound of the summation terms to be infinite. From the 
leading-order behaviour (12.101) . we know that the solution is singular at s = 0, and hence that 
X = 0 at s = 0. Solving (El with this boundary condition gives 


X = ± 


log(-2 


x/3) + 2Mtt\ s. 


where M is some integer. As in Type A solutions, the late-order behaviour for u n will take 
the form of an infinite sum of factorial-over-power ansatz terms, each associated with a unique 
choice of x- However, at late orders, this sum will be dominated by the four terms associated 
with M = 0 and M = — 1, which have the smallest absolute value. Hence, the four terms that 
contribute to the exponential behaviour of the solution are 

Xi = l°g(— 2 - \/3)s, X2 = -log(-2 - V3)s, 

X3 = log(-2 + \/3)s, X4 = -log(-2 + \/3)s. 


We observe that, unlike Type A, there are four sets of late-order behaviour to be considered in 
order to understand the Stokes structure of the associated solution. By solving (12.14a[) - (|2.14bl) . 
we again find that U' = 0. However, in this case, there are four distinct singulant expressions, 
and hence four constants associated with the corresponding late-order terms (denoted Ai - A4). 
We also find that, as before, k = —1/2. 

Hence, we find that the late-order terms may be expressed as 

Air(m/2-l/2) , A 2 r(m/2-l/2) 

Um ~ Vm ~ (log(—2 - v / 3)s) m / 2 - 1 /2 + (-log(—2- v^s)™/ 2 - 1 / 2 

A 3 r(m/2 - 1/2) | A 4 r(m/2 - 1 /2) 

+ (log(—2 + v^s)™/ 2 - 1 / 2 + (- log(—2 + v^s)™/ 2 - 1 / 2 ’ 

(4.2) 

where Aj is associated with Xi- 

The numerical process to determine Aj is more complicated than that demonstrated in Appendix 
El as we must now compute four values of Aj, rather than two. An example numerical calculation 
is illustrated in Appendix IA.2I 


4.2 Stokes structure 

Having calculated the singulants, we may determine the Stokes structure in a similar fashion to 
Section [3721 noting that Stokes lines follow curves on which Im(%) = 0 and Re(x) = 0, and anti- 
Stokes lines follow curves satisfying Re(x) = 0. The behaviour of the four different singulants 
is illustrated in Figure mi 
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Im(s) 



Im(s) 



Figure 4.1: Behaviour of the singulants for Type B solutions in the complex s-plane. The dashed lines indicate 
where the real part of the singulants is zero (and hence corresponds to an anti-Stokes line), and the 
solid lines indicates where the imaginary parts are zero. The jagged line represents a branch cut. 
Recalling that Stokes switching may only occur when Im(x) = 0 and Re(x) > 0, we see that \2 and 
X 4 switch across the Stokes line in the upper-half plane, while xi and X3 switch across the Stokes 
line in the lower-half plane. 
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From this figure, we see that the Stokes structure of the solution is substantially more com¬ 
plicated than in Type A solutions. Not only are there four different exponential contributions 
to consider, but there are two different sets of Stokes lines, each associated with two of the 
singulants. 

There is, however, also a more interesting difference. From Figure [TTTT1 we were able to conclude 
that neither exponential contribution could be present on the positive real axis in Type A 
solutions, as they would not be exponentially small. Hence, as the jump across Stokes lines is 
prescribed in (13.7D - (13.8D . the asymptotic behaviour is specified uniquely, with no free parameters. 
However, from Figure mi we see that this is no longer true for Type B solutions. While we can 
see that \2 and X 3 cannot be present on the positive real axis, as they would be exponentially 
large in the asymptotic limit, the same cannot be said about xi and X' 4 - As Re(xi, 4 ) > 0 on 
the positive real axis, the associated exponential contribution is exponentially small, and the 
asymptotic series is still valid when such contributions are present. Hence, the value of Six on 
the real axis may be freely specified, meaning that the solution contains two free parameters. 

Consequently, we determine that the exponential contributions active in each region of the plane 
are those shown in Figure B~2I ab We see that, assuming that dm 4 is nonzero on the real axis, 
the region in which the asymptotic expansion associated with (12.101) is valid is a sector between 
the anti-Stokes lines associated with xi ,2 and X' 3 , 4 ) in which the exponentials associated with 
X '2 and X 3 are inactive, and those associated with xi and X 4 are exponentially small. The angle 
between the edge of the sector and the real axis (denoted by 9 in Figure H~iH) is given by 


9 = tan 1 


^ log(2 + V3) ^ 


22.74°. 


However, we note that it is possible to choose (Six = 0 on the real axis. Figure l4~2T b1 illustrates 
this behaviour. In this case, the region in which the asymptotic expansion is valid is substan¬ 
tially larger than the general case shown in Figure H~2T aj: there is no exponential contribution 
associated with xi or X 4 to become large as the anti-Stokes lines are crossed. Instead, the angle 
9 is now given by 


9 = 7 r — tan 1 


/ log (2 + y^) \ 


157.26°. 


4.3 Stokes smoothing 

Using a nearly identical analysis to Section HOI we are able to show that for Type B solutions, 
the quantity switched across each Stokes line takes the form 

R n ~ SUe~ x/e , (4.3) 

where Rn is again the exponential contribution to the solution, and S is the Stokes switching 
parameter. Using matched asymptotic expansions, we find that 

K1+ „ 6n 

x' sinh(-x') ‘ 

We note that 

X' 1,2 sinh(—x'i )2 ) = V / 31og(—2 - V^), *3,4 sinh(-x' 3 x) = -log(-2 + V^). 

Hence, from (14.31) . the full exponential behaviour takes the form 

R n ~ »S 1 A 1 (—2 - V3)~ s/e + S 2 A 2 (—2 - V3) s/e 

+ S 3 A 3 (—2 + Vz)~ s,e + S 4 A 4 (—2 + Vz) s/e , 
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(a) General Stokes Structure 
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(b) Special Asymptotic Solution 


Figure 4.2: Behaviour of the singulants for Type B solutions in the complex s-plane. We illustrate the behaviour 
of the four exponential contributions (associated with the four singulants) individually, and also the 
composite behaviour when each is taken into account. The region in which the series expansion 
J2.6H provides a valid asymptotic approximation for u(s , e) and u(s, e) is illustrated in the composite 
expression. In each figure, the change in Stokes multiplier S across Stokes lines is illustrated. We 
recall that the jump across Stokes lines is specified in (I4.4D . Hence, as £ 2,3 must be zero on the real 
axis, their behaviour is completely specified. However, this is not true of <Si t 4 , which leads to two 
free parameters being present in the problem (specifically, the value of these Stokes multipliers on 
the real axis). Figure (a) illustrates the behaviour for general, non-zero choices of «Si ,4 on the real 
axis. Figure (b) illustrates a special case when <Si ; 4 are chosen to be zero on the real axis, giving a 
larger region of validity for the asymptotic solution. 
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where Si switches rapidly across the relevant Stokes line such that the difference is given by 


Pi, 2 ]t 


27r v / 3 


l°g(—2 - v'S) ’ 


[53,4]± 


—27t v / 3 

l°g(—2 + V3 )' 


(4.4) 


We also note for completion that S 2 and S 3 must equal zero on one side of the Stokes line, and 
hence are uniquely determined, while this condition is not required for 5i and £4. However, if 
both S+ and S£ are zero, we obtain the special Stokes structure illustrated in Figure POT bl. 

Consequently, the optimally-truncated asymptotic series for it(s, e), and hence v(s,e), is shown 
to take the form 


N 0 pt 

u(s, e) ~ e m/2 «m(s) + 5 1 A 1 (—2 - Vz)~ s/e + S 2 A 2 (-2 - \/3) s/e 

m =0 

+ 5 3 A 3 (-2 + Vl)~ s/t + S A A 4 (-2 + Vz) s/e (4.5) 

where the leading orders of u n are give in (12.10|) . and the late-order behaviour is given in 
m ■ This expression is an accurate asymptotic approximation up to exponentially-small terms. 
Unlike the Type A solution, however, this approximation of the Type B solution is only valid 
within a wedge of the complex plane illustrated in Figure 14.21 We also recall from Section 14.21 
that 1S2 and S 3 are uniquely determined, while <Si and <S 4 contain a Stokes-switching parameter 
that may be freely specified; this behaviour is not present in the asymptotic approximation of 
Type A solutions. 


5 Conclusions 

In this paper, we used methods of exponential asymptotics to compute and investigate two types 
of asymptotic solutions to the first discrete Painleve equation, associated with the leading order 
behaviours described in (12.91) and (12.101) . We identified the Stokes structure of the solution 
for each type, and used this information to deduce the sectors displaying asymptotically pole- 
free behaviour in the far-field limit. These asymptotic solutions were given as the sum of a 
truncated asymptotic power series and an exponentially-subdominant remainder term, shown 
in (13.9I) - (I3.10D for Type A solutions, and (14.51) for Type B solutions. 

In Section O we showed that the asymptotic solution associated with (12.91) is a valid approxima¬ 
tion in the far field of the entire complex plane, aside from a branch-cut along the negative real 
axis, and hence the solution is asymptotically pole-free. This solution has no exponentially-small 
behaviour on the positive real axis, although there is exponentially-small behaviour present in 
the left-half complex plane. This implies that this asymptotic behaviour is uniquely defined, 
with no free parameter present in the solution. 

The asymptotic solution associated with (12. 101) , discussed in Section 3J is generally asymptot¬ 
ically pole-free within a narrow wedge in the complex plane containing the positive real axis, 
seen in Figure !4~(?T al. Unlike the first asymptotic solution, however, there are two free param¬ 
eters in the asymptotic approximation, associated with Stokes switching. If these parameters 
are chosen appropriately, we discover asymptotic solutions with larger ranges of validity, such 
as that shown in Figure [j^ b). Outside of the asymptotically pole-free regions, the asymptotic 
series solution is no longer valid, and we instead see the elliptic function behaviour described 
by 0. In fact, the oscillatory behaviour that is present along anti-Stokes lines is a degenerate 
form of this behaviour. 

We note with interest that the two classes of asymptotically pole-free solutions to dPi, Types A 
and B, demonstrate power series-type behaviour in the far field, however they differ in that Type 
A solutions are uniquely specified with no free parameters, while Type B solutions contain free 
parameters associated with the Stokes multipliers in the solution. Similar asymptotic behaviour 
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is seen in the tritronquee and tronquee solutions to Pi, the former of which is uniquely specified, 
while the latter contains free parameters associated with Stokes switching. Consequently, the 
two classes of asymptotic behaviour described in this study may be regarded as corresponding 
to tritronquee and tronquee solutions of Pi respectively. 

The asymptotic solutions to dPi also provide insight into the behaviour of solution hierarchies 
to the equations described in m and CSD. The first of these describes solution hierarchies 
to the fourth Painleve equation. For example, Type B solutions (12.101) correspond to solution 
families which grow as w„(z) ~ ± \J—n/ 3 — z/Z + 0(n -1 / 2 ) in the large-n limit, with the 
parameters described in m also growing appropriately. These solution hierarchies are of 
particular mathematical interest due to their power series dependence on n in the asymptotic 
limit, which is particular to the asymptotic classes of dPi solutions described in the present 
study. 

Finally, we emphasize that this methodology may be applied to a wide range of difference 
equations, as it does not require the problem to be integrable, nor does it require that solution 
be expressed in terms of an integral equation. This is demonstrated in Appendix [Cl where we 
outline how this procedure could be applied to determine the Stokes switching behaviour of a 
non-integrable difference equation. Consequently, although we applied the method to determine 
information about the exponentially-small components of the solution to dPi, it could easily be 
adapted to analyze asymptotic solutions of broad classes of difference equations. 

One particular case of interest is the behaviour of orthonormal polynomials satisfying (II.dl) for 
a variety of different weight functions, such as the semi-classical weights described in [2]. In 
general, these polynomial families may be computed using non-integrable difference equations. 
However, as there is no aspect of the methodology which depends on the integrability of the 
equation, these difference equations are still amenable to the asymptotic methods described in 
the present study. 
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A Calculating the prefactor constants 


A.l Type A 


To determine the value of the constants A$, we return to the late-order expression given in (13.51) . 
We see that these expressions may be rewritten as 


Ai — A 2 (—l) m / 2-1 / 2 


M m (i7rs/2) m / 2 1 / 2 

r(m/2- 1/2) 


A 3 - A 4 (-l ) m/2_1/2 


t> m (i7rs/2) m / 2 1 / 2 

r(m/2 — 1/2) ' 
(A.l) 


By adding (and subtracting) successive terms appropriately in the limit to —> oo, we may easily 
determine Ai and A 2 by considering sufficiently high values of u m . For example, adding and 
subtracting appropriate even terms gives 


2Ai 

— lim 

'u 2 m (i7rs/2) m 1 / 2 

M2m+2(i7Ts/2) m+1 / 2 


m—> 00 

r(m- 1/2) 

r(m+l/2) 

2iA 2 

= lim 

'zi 2 m (i 7 rs/ 2) m " 1 / 2 

'« 2 m+ 2 (i 7 rs/ 2 ) m+1 / 2 ' 


ra—Kx) 

r(m — 1/2) 

r(m + 1/2) 


(A.2a) 
(A.2b) 


Similar results may be obtained for A 3 and A 4 by considering the behaviour of v n . 


As a specific example, we consider the case when a = — 1 , /3 = 1 and 7 = 0 . A numerical 
study similar to that performed by Joshi and Kitaev [2D], shown in Figure [Adi ah indicates that 
this behaviour tends to the appropriate asymptotic expression with the boundary conditions 
w 0 = 0.52040003 and wi = 0.55549107. We select the case when uo = \fs, giving u\ = 0 and 
112 = We may now evaluate (I2.11al) - (l2.11bl) in order to determine progressively higher 

terms in the series. In Figure lA.lf bl . we apply (IA.2al) to approximate Ai for various values of 
m, with the dashed line indicating the quantity that Ai converges to for sufficiently large to. 
Using the general relation in (I2.11al) - (l2.11b() . we calculate behaviour of u m up to to = 251, we 
find 


Ai « -0.08986 - 0.08986i, A 2 « -0.08986 + 0.08986L 

As 7 = 0, we have u m = —v m , we also conclude that A 3 = —A 3 and A 4 = —A 2 . This limit may 
be computed for different values of a and /3, although convergence can be somewhat slow for 
some choices of parameters. Finally, if we permit 7 to take nonzero values, we must consider the 
late-term behaviour of u m and v m independently, and solve for the four values of Aj separately. 


A.2 Type B 


When considering the more complicated form of the late-order terms for Type B solutions given 
in S3, we must solve for four different values of A.;. In this case, we write xi = |xi|e ie , and 
note the late-order terms may be rewritten as 


Ai - (—l) m/ 2 A 2 +e ie(m - 1 } A 3 - (-l) m / 2 e ie(m - 1 } A 4 


u m (log (-2 - \/3)s ) m/2 1/2 


(A.3) 


T(to/2- 1/2) 

While this expression is not as simple as (ED. we may still use it to determine the values of 
Ai by taking four distinct values of m and solving the resultant system of equations. As these 
values of to become large, the computed values for Aj will converge to the appropriate values. 
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(a) w n 
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Ai ss -0.08986 - 0.08986i 



• Re(Ai) 
o Im(Ai) 


25 
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(b) Approximation for Ai 


100 


Figure A.l: Figure (a) illustrates the behaviour of the solution to dPi with a = — 1, (3 = 1 and 7 = 0 . The 
boundary conditions are chosen such that wq = 0 and w\ = 1. The values of w n are represented 
as black circles, and the asymptotic solution w n = d zy/—an, associated with Type A solutions, is 
represented as a solid curve. It is apparent that the behaviour of the discrete equation tends to 
the asymptotic expression in question for large n. Figure (b) illustrates the approximation for Ai 
obtained by computing the expression given in (IA.2al) . As m increases, the approximation tends to 
the dashed curve. 


We consider a numerical example, choosing a = 3, /3 = 1 and 7 = 0 (ensuring that u m = v rn . 
and hence eliminating the need to consider v m separately). We begin with the leading-order 
solution zto = y/s and compute the behaviour of u m up to m = 500 and use four consecutive 
even values of m in (IA.3I) to generate the system of equations which is subsequently solved for 
The computed values of Ai and A4 are shown in Figure lA~2l for a range of values of m , 
and we see that they are converging to a constant value. By taking a system of equations using 
sufficiently large computed values of zz m , we find that 

Ai « 0.04666 + 0.02669i, A 2 « 0.11212 - 0.06491i, 

A 3 w 0.04666 - 0.02669i, A 4 fa 0.11212 + 0.06491L 

We have therefore determined the form of the late-order terms u m , and consequently v m . Again, 
if we considered an example with 7 / 0, we would find that u m ^ v m , and the values of A ? ; 
associated with v m would need to be computed separately. 


B Stokes smoothing 

To perform an exponential asymptotic analysis on the remainder terms in (13.61) , we must ensure 
that the asymptotic series is truncated optimally. To find the optimal truncation point, we follow 
the commonly-used strategy described by m , in which the series is truncated at its smallest 
term (and hence where consecutive terms have the same size). 

By solving 

e ( iv + 1 )/ 2 Ujv - +1 

e N u N ~ ’ 

in the limit e —> 0 and N —> 00, we find that N ~ 2|y|/e. Hence, we will set 

7V=^M+2w, (B.l) 
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(b) Approximation for A 4 


Figure A.2: This figure illustrates the computed approximations associated with Type B solutions having a = 3, 
(3 = 1 and 7 = 0 for (a) Ai and (b) A 4 , obtained by solving the system generated from fOt using 
four consecutive even values of m, beginning with the value indicated on the horizontal axis. As m 
increases, the approximations tend to the dashed curve. 


where 0<w<l/2is chosen so that N takes integer value; the choice of the two in 2 u is one of 
algebraic convenience. 

Applying this series expression and using (12.11al) - (|2.11bl) to eliminate terms, recalling that uq = 
-vo, gives 


2uqSn + 2e 1 /“uiSjv + 2eu 2 SN + 2ug ^ ' 


L ^ 2J 


7=1 


m 


LJV/2J-1 2 • (2 j) LJV/2J-1 2 ■ (2 j) 


LAT/2J £ 2j^j(2j) 

2vqRn + 2e 1 ^ 2 v\R n + 2ev 2 RN + 2vo ^ A 

3 =1 


(2.7)! 


UV/2J-1 2j p (2j) LA/2j-l 2j (2 j) 

2^/V E -$r + 2 ™ 2 E N 


7=1 


7=1 


m 


+ 2e N / 2 voiiN + ... (B.2a) 


e 1 ^ 2 7 Rn + 2e N / 2 uoi>N + ■ ■ ■ (B.2b) 


where the omitted terms will be smaller in the limit e -> 0 than those which have been retained. 

Noting that Rn and Sn have decoupled, we will use (IB.2bl) to determine the behaviour of Rn- 
We apply a WKB ansatz to the homogenous version of (IB.2bl) by setting Rn = /(s)e s ^ s ^ e . The 
solution away from the Stokes line (where the right-hand side not have any effect) is therefore 
given by the homogenous version of (IB.2bl) . and therefore to leading order as e —> 0 


By comparing this to the singulant equation JO, we see that we can write g(s) = %, where \ 
represents the value of the singulant under consideration; here, the analysis will be performed 
in the general case, with the particular choice of \ made afterwards. Continuing to the next 
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order in e, we will find that f(s) is given by a constant. Hence, away from the Stokes line, the 
exponentially-small remainder takes the form Rn ~ Ce~ x ^ e , where C is a constant, which we 
will set to be the associated value for A. 

Knowing that the behaviour will switch rapidly in the neighbourhood of the Stokes line, we set 

R n = SAe-*/ e , 


where S is the Stokes switching parameter that varies rapidly as the Stokes line is crossed. The 
multiplier A is chosen for subsequent algebraic convenience. Applying this ansatz to (IB.2al) 
gives, after some cancellation, 


L 7 J Lv'l2j-l 

2 «' S ' A E 

As N —> oo, the series may be taken to be infinite, introducing only exponential error. However, 
we have already discarded terms containing algebraic error, so this is not a concern. We therefore 
write 

<S'Asinh(—x') ~ e ^/2- 1 e x/e Ujv _ 

Noting the form of un given in (12.121) . we write 


d5 ^ jy/2-1 y/e W/2 + k) 
ds sinh(— x')x N ^ 2+k 


(B.3) 


Noting the form of the optimal truncation m, we set x = re 16 and restrict ourselves to curves 
with r fixed. This gives 

d ix'e -1 ® d 

ds r dO ’ 

while (IB. II) becomes N = 2r/e + 2 u. Hence (IB. 31) becomes 

— „ rc i 8 c r,e+UJ ~ 1 r(r/e + uj + k) 
dd x' sinh(—x , )(re I0 ) r / e+a, + fc 

From (13.41) . we note that x' sinh(—x') = 7r/2 for both xi and X‘ 2 - We also recall that k = —1/2. 
Hence, applying Stirling’s formula and simplifying eventually gives 

- i/^exp(-(e ie -l-i0)-i0(A : +w-l)). 

This is exponentially small, except in the neighbourhood of 9 = 0, which corresponds to the 
Stokes line (that is, x being purely real and positive). To determine the local behaviour, we 
apply the scaling 6 = e 1 / 2 #, which gives 


Hence, by integration, we find 


5 




rO/VT 


- fi2 / 2 ds + C 


(B.4) 


where C is an arbitrary constant. Importantly, as the Stokes line is crossed, 9 varies rapidly 
from —oo to oo. Hence, away from the Stokes line (0 = 0), we find that the jump in the Stokes 
multiplier A as the Stokes line is crossed is given by 


[5] 



= 4, 
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and hence that the exponential contribution switched across the Stokes line Rjy changes by 

[R n ]± ~4Ae- x/e , 

as s moves from negative to positive. Recalling that there are exponential contributions from 
the two singulants given in (USD, we find a combined expression 

Rn ~ SiAi e -^/ 26 + S 2 A 2 e^ 2 \ (B.5) 

where 5i changes rapidly from four to zero when crossing from the third to the fourth quadrant 
in the smooth manner described by (ES, and S 2 changes from negative four to zero when 
crossing from the second to the first quadrant in the same way. A similar analysis shows that 

R n ~ SiA 3 e -' ins/2t + S 2 A 4 e™ s/2e . (B. 6 ) 


C A non-integrable equation 


To demonstrate the utility of this method, we consider a non-integrable variant of dPi, and 
briefly outline the steps required to determine the late-order term behaviour, and hence the 
Stokes structure of the solution. We begin with the equation presented in (EH), and adjust the 
final term, giving 


w n+ i + w n + w n -1 


an + /? 

-h7®„. 

W n 


(C.l) 


Representing this as a system of equations in y n and z n . defined as in EH- We again using the 
method of multiple scales by considering the slow variation in terms of s == en, and apply the 
transformation given in El- Hence the rescaled equations become (after Taylor expansion in 

s) 


* 0 +) + E 


3 =0 


^ 0 +) + E 


j=o 


2e 2j v( 2j )(s, e ) 

m 

2 e 2 ^u^(s, e ) 

m 


u(s, e) = as + e/3 + 7 u(s, e) 2 , 
v(s, e) = as + e/3 + jv(s, e) 2 . 


We expand the solutions u{s,e) and v(s,e) as a formal series in powers of e, rather than half¬ 
integer powers as in section [ 2 j as there is no e 1 / 2 present in the rescaled equations. We denote 
the coefficients multiplying e m as u m and v m . As in the analysis of dPi, there are two sets of 
solutions, associated with the leading-order behaviours 


uo = -Vo = ±J- 


1 + 7 


u 0 = v 0 = ± 


3-7 


which we again refer to as Type A and Type B solutions respectively. We find that the general 


relation for higher- 

-order terms becomes 
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(C.3a) 
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(C.3b) 


3 =0 

3 =0 

r= 0 

3 =0 



We again define a late-order ansatz in the limit n —> 00 , however due to the different form of 
the asymptotic series, we instead select the ansatz to have the form 


U(s)T(m + k) V(s)T(m + k) 


(C-4) 
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Applying this ansatz to IIC.8all - IIC.3bll and matching leading orders of m as m —> oo gives new 
singulant and prefactor equations. The singulant equations for Type A and Type B solutions 
respectively are now given by 

Type A: cosh (—= 7, Type B: 2 + cosh(— x') = 7, 

with the condition that x(0) = 0 in both cases. By considering higher orders in m as m —> oo, 
we may obtain new prefactor equations, and perform a Stokes switching analysis in a similar 
fashion to the analysis performed on dPi. 

We therefore see that, although the exponential asymptotic method described in this study may 
be used to study the Stokes-switching behaviour of integrable equations such as the discrete 
Painleve equations, it does not depend on the integrability of the problem, and hence may be 
applied to a much more broad class of equations without significant changes to the process. 
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